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Abstract 

We investigate the effect of exciton coupling on the optical absorption spectrum of polymer 
molecules under conditions of strong inhomogeneous broadening. We demonstrate that the de- 
pendence of the maximum in the rescaled absorption spectrum on the number of monomers is 
determined by the average monomer excitation energies and their resonant coupling and insensi- 
tive to the inhomogeneous broadening. Thus the absorption spectrum can be used to determine 
optical interactions between monomers. The results are applied to the absorption spectra of poly- 
A poly-T DNA hairpins and used to interpret the dependence of the absorption spectrum on the 
number of monomers. We also discuss exciton localization in these hairpins. 

PACS numbers: 7080.Le, 72.20.Ee, 72.25.-b, 87.14.Gg 
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1. The optical absorption spectra of dissolved molecules are strongly affected by in- 
teraction with the environment. The positions and orientations of solvent molecules with 
respect to each molecule differ, and leads to different shifts in their optical excitation ener- 
gies. These shifts can be characterized by a distribution with characteristic width W that is 
usually called inhomogeneous broadening. In many important systems, particularly includ- 
ing biological molecules in their native environment involving a highly polar solvent (water), 
this broadening is very large W > O.leV and it can smear out various interesting features 
of the spectrum including its sensitivity to molecular conformation and geometry. 

However there exist experimental indications that despite the large inhomogeneous broad- 
ening the absorption spectrum contains important information about molecular properties. 
Recent measurements of optical absorption spectra of DNA hairpins made of identical AT 
base pairs demonstrate a remarkable change in absorption maxima with the number of 
monomers []| . We believe that this effect is due to sensitivity of the electronic excited state 
to the number of monomers. The excited state (exciton) can be delocalized between all n 
identical monomers and the resulting exciton energy (absorption energy) is sensitive to the 
size of the overall wavefunction. The exciton coupling Vij of different monomers is caused 
by their resonant interactions, and it reaches its maximum value Vq for adjacent monomers. 
One can therefore expect that the change in the absorption spectrum with n is determined 
by Vq. Since the multipole interactions that determine Vq are very sensitive to the distance 
between monomers and to their relative orientation, the change in absorption spectrum with 
n can be used to study DNA conformation and structure Q]. Thus time-resolved optical mea- 
surements of DNA participating in biological or chemical processes like protein attachment 
can provide important information about the real time kinetics of such process and DNA 
structural and conformational changes during the process. In addition, the understanding of 
exciton coupling is important for DNA applications as a naturally available optical material. 

Below we investigate the inhomogeneously broadened spectrum of a polymer molecule 
made of n identical monomers assuming that inhomogeneous broadening W exceeds both 
the exciton coupling Vq and the variation A of the average excitation energies of the different 
monomers with respect to each other caused by their different positions (e. g. monomers 
at the edges compared to monomers in the middle of the polymer). Under these conditions 
the absorption spectrum a(E) has a single maximum corresponding to electronic excitation. 
The position of the maximum is the important measurable characteristic of the spectrum 



2 



which can be easily extracted from the experimental data. It turns out that the more 
informative characteristic is the maximum in the rescaled absorption /3(E) = a(E)/E. We 
investigate the dependence of the maximum in (3(E) on the number n of monomers, showing 
that in the limit of strong inhomogeneous broadening this dependence is determined by the 
exciton coupling Vo and the monomer energy variation A, while it is insensitive to the 
inhomogeneous broadening. If fluctuations in the monomer energies are independent and 
obey Gaussian statistics then the change in the maximum position is determined by the 
simple analytical expression 

SE_ = (!) 
n 

where 9 is the angle between the transition dipole moments of adjacent monomers. As we 
show below this analytical result with the parameter taken to be Vq ~ 0.042eV, as obtained 
from semiempirical calculations and the DNA twisting angle 9 = 36°, agrees with the results 
of recent experiments (Ref. l| and Fig. [T]). The localization of excitons in DNA will be 
discussed in the light of these findings. 

2. Optical excitation in a chain of n monomers can be described by a tight binding Hamil- 
tonian coupled to the environment. Differences in the environment for different monomers 
lead to fluctuations in the monomer excitation energies fa 

n 

H = J2v ijC t Cj + J2^i4^ (2) 

i^j t=l 

Here q, cf are operators for creation and annihilation of an exciton in a site (monomer) 
i. The excitation energies fa are taken to have random values that are characterized by 
a distribution function P(fa, 4>2, ■■<pn) and they are directly coupled to the local exciton 
population n, = cfci. All monomers % = 1,2, ...n possess transition dipole moments Hi 
having identical absolute values /i - Since the photon absorption process is very fast one 
can the treat energies fa as static, because their characteristic time of change is defined 
by a slow atomic motion. For the same reason we only consider vertical excitations, as 
the environment cannot relax to the new configuration during the short absorption event. 
We ignore fluctuations of exciton couplings due to environmental fluctuations as this 
interaction has a power law dependence on interatomic distances that results from multipole 
interactions. Since fluctuations in the interatomic distances are much smaller than the 
distances themselves, the power law dependent interactions should fluctuate weakly. 



The position of absorption maximum depends on the distribution function P(fa, 02, ■•4>n)- 
Below we assume that this energy distribution function can be approximated by a Gaussian 
distribution. This is a reasonable assumption when fluctuations of energy are smaller than 
the energy itself which justifies the use of a second order expansion of the free energy for 
the probability distribution P [4]. On the other hand the width of the distribution (inho- 
mogeneous broadening W ~ 0.2eV) is much larger than the natural linewidth hT < 10~ 3 eV, 
which allows us to ignore deviations from the Gaussian lineshape due to Lorentzian tails in 
the natural line. Of course calculations can be easily extended to any other distribution P. 

If the distribution P is a Gaussian distribution then all other distributions for the reduced 
number of energies obtained by the integration of P with respect to the remaining energies 
are also Gaussian distributions. For our study we will need distributions either for a single 
energy or for two energies. We introduce the single energy distribution as 

m)-^^e^(^-j/), (3) 

where fa is the average excitation energy of the monomer i and Wi is its dispersion (inho- 
mogeneous broadening). The distribution function for two energies can be expressed as 

2W 2 (1-k 2 .) 2W 2 (1-k 2 .) Wi W,-(1-k?.) 

Piitt, <P') = t— > (4) 

where the parameter Kij describes correlations of the excitations of monomers % and j. Note 
that < (fa — fa)(fa — 4>j) >= K ijWiWj. We further assume that all average energies fa 
are close to some typical energy O an d their deviations from it are small compared to the 
inhomogeneous broadenings Wi. The same is assumed for exciton coupling so for further 
considerations we set 

Va,\fa-fa\<^Wi. (5) 

Here we have assumed that the energy gap to the next excited state exceeds the inhomoge- 
neous broadening so we can neglect it. 

3. One can describe the optical absorption of the ensemble of polymers in terms of exciton 
eigenstates ip a >,a = 1, ...n, their energies E a , transition dipole moments \i a and oscillator 
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strengths f a (Eq. <^j) as 

n 

a(E) oc < fAK ~ E) >, 

0=1 

_ 2fx 2 a mE a 

Ja 3(eh) 2 ' K) 
where e and m are the charge and mass of electron. 

Consider the exciton coupling V as a perturbation. In the zeroth order approximation in 
V the eigenstates a coincide with individual monomer i excitations characterized by energies 
Ei = <pi and all the transition dipole moments are identical, /ij = fiQ, so we get 

> <xj2 d<f>p t (mE-<P) = J2 p ^- ( ? ) 
»=i i=i 



To first order in V, which we are interested in, there are no corrections to the exciton 
ene rgy , because the oMagcal faction V ^ ie S eco, ld o r de r colons o, lly & 
while there is a correction to the transition dipole moment associated with the change in 
exciton wavefunction \i >= \i > + \ j > This correction can be expressed 

as /Ij, — fj, i + J2j J Lil- l jVij/(,Ej — E^. Accordingly the change dfi in the oscillator strength 
can be expressed as 5 fa/ fa = 2 Y^j^i cos (%)^ij / '{Ej ~ Ei), where is the angle between the 
transition dipole moments of the monomers i and j. Finally the rescaled absorbance can be 
expressed as 



(3(E) = ^l^j2 p i( E ) + 2 E v * cos ( 

1=1 ijtj 



xRv - L # ^V- (8) 

The factor E in the definition of a(E) makes the absorption maximum sensitive to inho- 
mogeneous broadening. To study lineshape effects it is more convenient to use the maximum 
in the rescaled absorbance (3(E) = a(E)/E, which represents the spectral density. Consider 
the energy corresponding to the maximum in j3{E). In the zeroth order approximation one 
can find this energy using only the first term in Eq. (jHJ). Under the conditions specified 
by Eq. (jHJ) one can expand all exponents in this term as e x ~ 1 + x and then set the 
derivative of the expression with respect to energy E to be equal to zero. Then we get 
J2ti(E {0) ~ <Pi)/Wf = 0. This yields 

=< 4> > w , (9) 



where < A > w = (£™ =1 A/Wf) / £? =1 

One can approximate the shift of the maximum the due to the second term in Eq. (jSJ) as 
« [dfi x \E^/dE^/[<PpW(E { ®)/d&\ t where Z^ ' 1 ) represents the first (or second) terms 
in Eq. (jSJ). Both first and second derivatives can be evaluated at E^ /W ~ 0. They can be 
evaluated as d^(E^)/dE « (2/V2tt) Ei# ^j7(W?Wj) where ^ = ^ cos(0 i>3 -)/(l + «y) 
and d 2 p i0 \E^)/dE 2 = (l/v 7 ^) ELi 1 / 1 ^ 3 - Accordingly the position of the maximum 
^ roai of the function (3(E) can be expressed as 



^mai =< >W H 1 7w3 " ^ iU J 



This result does not depend on the absolute value of the inhomogeneous broadening. 

In the simplified approximation one can consider only nearest neighbor interactions Vq. 
Thus we assume Wi = W n = W ex , Wi = Wi nt , i = 2,3, ...n — 2 and 0i = (p n = </>o + A, 
(f>i = 4>a, i = 2, 3, ...n — 1. Then Eq. ( JTU1) can be rewritten as 



^max = 00 + 4K + 
2A _A\r( 1 _L _2 1 1 

n-2 j 2_ 



(11) 



where Vo = V \ c "^ , is the correlation coefficient for adjacent monomers. This result is 
different from the maximum shift obtained using the convolution by Gaussians method 6| 
because the latter method ignores the relative energy fluctuations of monomers. 

4. It is interesting to compare our predictions with experimental data. Here we consider 
the measurements of light absorption by poly-A - poly-T DNA hairpins endcapped by stil- 
bene jl| and containing n AT base pairs. To analyze the data we extracted positions of 
spectral maxima in /3(E) for n = 1, 2, 3, 4, 5, 6, 8, 12. We used the second maximum from 
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the low energy side because the first maximum is due to stilbene absorption 

Since all hairpins contained stilbene groups which also contributed to absorption the 
definition of maximum gets more accurate at large n. Therefore we fit experimental data 
using Eq. (TTTT) in the form E max = E OQ —a/(bn+c) for n = k, k+1, ...12 choosing k = 1,2, 3, ... 
The results of the data fit become almost insensitive to our choice of k for k > 3 so we choose 
k = 3 as a lower constraint. For smaller k fitting parameters can change by 50%. For any k 
the best data fit is obtained setting c = (W ex = W in ) in Eq. (TTTT) which suggests that the 
fluctuations of base pair energies are identical for DNA bases located at hairpin edges and 



inside the sequence. The reasonably good data fit is obtained setting (see Fig. [T]) 

Emax(n) = 4.5242 - 0.135/n. (12) 

One can assume that the energy difference A (Eq. (jlip ) of end and middle hairpin 
AT sites is negligible, e. g. due to pronounced screening effect of the static interactions 
leading to the "diagonal" energy difference A. For dynamic "off-diagonal" interaction Vq the 
screening is much smaller because it has to be taken at high exciton frequency u > 10 15 s _1 . 
Then Eq. ([IT]) can be rewritten as (cf. Eq. ([1])) 

W cos(fl) 4V cos(fl) 

E ™=^ + 1 + KQ - n (l + K )- (13) 

Here we have assumed 9 = 36°. Then the agreement of theory and experiment requires Eq. 

Vo/(l + «o) ~ 0.0417eV. (14) 

To verify this expectation we calculated the exciton coupling strength Vq using the semiem- 
pirical program ZINDO (part of the Gaussian 03 package .8]). The input geometry was 
obtained using 3DNA software (7J with several AT base pairs placed at the neighboring dis- 
tance 3.4A with the twisting angle 36°. The ground state input geometry for each AT pair 
has been optimized using density functional theory calculations based on B3LYP/6-31G** 
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■ Experimental data [1] 
Optimum data fit using Eg. (11) 
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FIG. 1: Theoretical fit of experimental data [l| for the dependence of maximum in the scaled 
absorption intensity (3(E) on the number of monomers for DNA hairpins (AT) n . 
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. We studied the first "bright" excited state of the AT pair [9|. The energy of this state 
is found to be 4.2eV. This energy is smaller than the experimental peak in the absorption 
spectrum at 4.5eV by 0.3eV which is a typical error for ZINDO 10|. The exciton coupling Vo 
has been defined as half of the splitting of this bright excited energy level for two AT pairs 
similarly to Ref. The extracted value Vo ~ 0.042eV is very close to the experimental 

prediction Eq. (jT4l) if we ignore correlations between the excitation energies of adjacent base 
pairs (/«o <C 1). Since these correlations are due to the long-range electrostatic interactions, 
they should be weaker for neutral excitations than for charges where k « 1/2 12j. Our 
results agree qualitatively with the analysis of excitons in the single strand DNA made us- 
ing A bases that were reported in Refs. 0, Q. Excitonic coupling has been estimated as 
0.053eV in that system. 

It is important that the exciton coupling Vq is positive leading to an increase of the 
absorption maximum energy E max with increasing number of monomers, i. e. a blue shift 
of the absorption spectrum in agreement with the experiment . This is because of the 

dipolar nature of exciton coupling. The oscillator strength is larger for quantum states with 
larger transition dipole moments. These states have higher energy because the interaction 
of parallel AT pair dipole moments is repulsive in DNA geometry. The red shift obtained in 
Re, Q is probably the consequence of self-interaction errors in the TDDFT calculations 
that were used. 

5. We have investigated the inhomogeneously broadened absorption spectrum of a 
molecule made of n identical monomers. The shift of the maximum of the rescaled ab- 
sorption /3(E) = a(E)/E with n is determined entirely by the average monomer energies 
and their exciton couplings. It is insensitive to the inhomogeneous broadening and can be 
expressed by a simple algebraic form Eq. (fTTl) . Thus exciton coupling can be characterized 
using inhomogeneously broadened absorption spectra. 

A related but challenging problem is whether the exciton quantum state is localized within 

nn 

a single AT base pair or it is extended through several base pairs [14], H5j . This localization 
is defined by the ratio of exciton coupling strength Vo and inhomogeneous broadening W. 
If Vq/W <C 1 this should be the regime of strong localization [4]. We can use our estimate 
Vn ~ 0.04eV and compare it to the width of the DNA absorption maximum W ~ 0.2eV 
jj. Then Vo < W so the exciton wavefunction is probably localized within almost a single 
site. Note that a recent study of this effect for single stranded DNA considered the behavior 
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of a large number of eigenstates, and came to the opposite conclusion; however this is not 
necessarily in contradiction with the present analysis, which only considered a single excited 
state More excited states need to be considered to verify our expectations. 
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